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ABSTRACT 

We report on new and archival X-ray and near-infrared (near-IR) observations of the anomalous 
X-ray pulsar IE 1048.1—5937 performed between 2001-2007 with the Rossi X-ray Timing Explorer 
(RXTE), the Chandra X-ray Observatory (CXO), the Swift Gamma-ray Burst Explorer, the Hub- 
ble Space Telescope (HST), and the Very Large Telescope. During its ^--^2001-2004 active period, 
IE 1048.1—5937 exhibited two large, long-term X-ray pulsed-flux flares as well as short bursts, and 
large (>10x) torque changes. Monitoring with RXTE revealed that the source entered a phase of 
timing stability in 2004; at the same time, a series of four simultaneous observations with CXO and 
HST in 2006 showed that its X-ray flux and spectrum and near-IR flux, all variable prior to 2005, 
stabilized. Specifically, we find the 2006 X-ray spectrum to be consistent with a two-component 
blackbody plus power law, with average kT = 0.52 keV and power-law index F = 2.8 at a mean flux 
level in the 2—10 keV range of ~6.5xl0~^^ erg cm~^ s~^. The near-IR flux, when detected by HST 
{H ^ 22.7 mag) and VLT {Ks ~ 21.0 mag), was considerably fainter than previously measured. 
Recently, in 2007 March, this newfound quiescence was interrupted by a sudden flux enhancement. 
X-ray spectral changes and a pulse morphology change, simultaneous with a large spin-up glitch and 
near-IR enhancement. Our RXTE observations revealed a sudden pulsed flux increase by a factor 
of ^-^3 in the 2—10 keV band. In observations with CXO and Swift, we found that the total X-ray 
flux increased much more than the pulsed flux, reaching a peak value of >7 times the quiescent value 
(2—10 keV). With these recent data, we find a strong anti-correlation between X-ray fiux and pulsed 
fraction. In addition, we find a correlation between X-ray spectral hardness and flux. Simultaneously 
with the radiative and timing changes, we observed a significant X-ray pulse morphology change such 
that the profile went from nearly sinusoidal to having multiple peaks. We compare these remarkable 
events with other AXP outbursts and discuss implications in the context of the magnetar model and 
other models of AXP emission. 

Subject headings: pulsars: general — pulsars: individual (IE 1048.1—5937) — stars: neutron — stars: 
pulsars 



1. INTRODUCTION 

The small handful of unusual objects known as anoma- 
lous X-ray pulsars (AXPs) are characterized by their 
slow spin periods (5-12 s), large inferred magnetic fields 
(~10^^ G), and persistent X-ray luminosities in great ex- 
cess of available spin-down power. These young, isolated 
pulsars have properties common to another class of ob- 
jects, the soft gamma repeater s (SGRs). For recent re - 
views, seelmods & Thompspjol ([2006) and lKaspil (|2007D . 
Potentially the most intriguing and most revealing prop- 
erty shared by AXPs and SGRs is their highly volatile 
nature: both are known to exhibit occasional and sudden 
dramatic flux and spin variability in the form of X-ray 
bursts, flares, and glitches. 

The magnetar model (iDuncan fc Thompso^ 119921 : 
[Thompson fc Duncan|[l995HThompson fc Duncanlll996D 

^ Department of Physics, Rutherford Physics Building, McGill 
University 3600 University Street, Montreal, QC, H3A 2T8, 
Canada 

^ NASA Goddard Space Flight Center, Astrophysics Science Di- 
vision, Code 662, Greenbelt, MD, 20771, USA 

^ NPP Fellow; Oak Ridge Associated Universities, Building 
SC-200, 1299 Bethel Valley Road, Oak Ridge, TN, 37830, USA 

Canada Research Chair; Lorne Trottier Chair; R. Howard Web- 
ster Fellow of CIFAR 

^ Dynetics, Inc., 1000 Explorer Boulevard, Huntsville, AL, 
35806, USA 

^ NSSTC, 320 Sparkman Drive, Huntsville, AL, 35805, USA 



was instrumental in explaining both the "anomalous" 
X-ray emission, as well the episodes of burst activ- 
ity. It identified the sources as highly magnetized iso- 
lated neutron stars, with field strengths on the order 
of lO^'' — 10^^ G. The magnetar model uniquely, and 
correctly, predicted that SGR-like bursts would be ob- 
served from AXPs. A major proposed effect of this 
enormous magnetization is the "twisting" of the magne- 
tosphere, with associated magnetospheric currents and 
heating of the crust at the twisted field-line footpoints 
(jThompson et al.]l2002l ). This gives rise to the unusual 
X-ray spectrum, which has generally been empirically 
well characterized by a two-component model consisting 
of a power law plus blackbody. Bursts of high-energy 
emission occur presumably when the crust succumbs to 
magnetic stresses and deforms, leading to a rapid re- 
arrangement of the external ma gnetic field. Recently, 
iBeloborodov fc Thompson! (|2007f ) proposed the existence 
of a plasma corona contained within the closed magneto- 
sphere to explain the broad band spectrum of magnetars 
that extends f rom the infrared (I R) to hard X-rays be- 
yond 100 keV (|Kuiper et al.ll2006( ). 

IE 1048.1—5937 has had an unusual history, even 
by AXP standards^. Prior to 2002, monitoring of 

^ A summary of its properties can be 
found in the online SGR/AXP Catalog 
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the pulsed flux with the Rossi X-ray Timing Explorer 
{RXTE) showed that its spin down was so unstable 
that phase coherence could b e maintained for peri- 
ods of only months at a time (|Kaspi et al.l |2001D , be- 
haviour which ec hoed earlier reports of P, flux and spec- 
tral variability (| Poster broek et all Il998f ). Two small 
bursts, the first ever found in an AX P, were dis- 
covered in 2001 from IE 1048.1-5937 (jGavriil et all 
[2002) , clinching the suspected association between AXPs 
and SGRs. Since 2002, tri-weekly RXTE observa- 
tions of IE 1048.1—5937 allowed the changing pulsed 
flux and spin-down rat e to be monito r ed on rela- 
tively short timescales. iGavriil fc Kaspil ()2004 here- 
after GK04) reported order-of-magnitude torque vari- 
ability on timescales of weeks to months that only 
marginally correlated with large luminosity variations. 
The unusually slow and long-lived pulsed flux "flares" 
(not to be co nfused with the SGR giant flares; see 



IWoods k. Thom pson 2006) that began in late 2001 and 
lasted into 2004 are in some ways unlike activity seen 
in any other AXP or SGR thus far. The 2-10 keV 
spectrum and pulsed fraction are also variable on long 
tim escales, the latte r being anti-correlated with luminos- 
ity (jMereghetti et a l. 2004; Ticngo et aL 2005). At near- 
infrared (near-IR) wavelengths, the flux varied dramat- 
ically during 20 02-2003 (.Wang fc Chakrabartvl l2002l : 
llsrael et"al1l2002l ): unfortunately, the sparsity of near-IR 
observations prohibits measurements of the variability 
timescale, althoug h a near-IR/X-ray anti-corr elation has 
been suggested (jPurant fc van Kerkwiikl l2005). Further 
activity in the form of anothe r burst and sm all pulsed- 
flux increase was seen in 2004 (IGavriil et al.l [2006'l. 

Here, we report on a series of simultaneous CXO and 
HST observations obtained during the course of 2006, as 
well as on archival CXO, XMM and VLT observations, in 
§ [2l We also present RXTE timing results from the past 
three years. In 2007 March, IE 1048.1-5937 unexpect- 
edly entered yet another new phase of activity, and we 
report on Target-of-Opportunity observations with CXO 
and Swift taken after the 2007 March event. In § 13.11 
we discuss the 2004-2006 period of radiative and spin 
quiescence in the context of both magentar and accre- 
tion models, and in § 13.21 we describe and interpret the 
events following the 2007 March flare. 

2. OBSERVATIONS, ANALYSIS AND RESULTS 

As part of a long-term project, we regularly monitor 
AXPs using the Rossi X-ray Timing Explorer (RXTE). 
IE 1048.1—5937 is by far the most frequently observed 
AXP due to its relatively poor rotational stability. The 
RXTE observations (see § 12. ip are crucial in measur- 
ing the spin evolution of the source and its pulsed 
flux. However, the high background and large unim- 
aged field-of-view of RXTE make spectral measurements 
of IE 1048.1-5937 diflrcult and total flux (and con- 
sequentially pulsed fraction) measurements impossible. 
To make these measurements, we obtained five Chan- 
dra X-ray Observatory observations {CXO; see § 12. 2p 
approximately equispaced throughout 2006. In order to 
probe the origin of the near-IR emission of this source, 
we observed with the Hubble Space Telescope {HST; see 
simultaneously with CXO. The motivation for si- 
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Fig. 1.— The long-term evolution of IE 1048.1-5937's pulsed 
properties. The three vertical dashed lines indicate the approxi- 
mate beginnings of the flux flares on (from left to right) 2001 Oc- 
tober 26, 2002 April 6, and 2007 March 21, respectively. Horizon- 
tal dotted lines represent the value of the plotted parameter prior 
to all observed flares, (a) Spin frequency as observed by RXTE. 
Small points show the measured frequencies, while intervals over 
which phase coherence has been maintained are shown as thick 
lines (GK04; Dib et al., in preparation). (6) Frequency derivative 
as a function of time. While the frequency was varying in 2002-2004 
and a phase-coherent timing solution could not be found, u was de- 
termined in short intervals b y calcul ating the slope of three consec- 
utive u measurements (after IGK04I) . Since 2004, phase coherence 
has been maintained and u reflects the standard technique (Dib 
et al., in preparation), (c) 2—10 keV rms pulsed flux as observed 
by RXTE. (d) Simulated total 2—10 keV unabsorbed flux, shown 
as green points. Total flux is estimated from the RXTE pulsed 
flux and the power-law correlation between pulsed fractio n and 
measured total flux, shown as red points, as described in § 13.2.21 
The simulated fluxes (green) have been scaled to match the ac- 
tual measured total fluxes (red), (e) 2—1 keV rms pulsed fra ction 
calculated using the method described in IWoods et al.l ()2004l ). 

multancous CXO and HST observations was to search 
for correlated variability in the different bands, as that 
provides insight into the physical mechanisms generat- 
ing the radiation. We also examined archival ifg-band 
observations with the Very Large Telescope {VLT; see 
§ 12. 4p , which contribute to our monitoring of its near-IR 
brightness. Following the 2007 March glitch detected in 
our RXTE moni toring observations of IE 1048.1-5937 
(jPib et al.ll2007bf ). we initiated three CXO and two Swift 
Gamma-Ray Burst Explorer {Swift; see ? l2.2p Target-of- 
Opportunity (ToO) observations to follow the X-ray flux 
and spectrum of the source. In order to characterize the 
long-term evolution of IE 1048.1—5937 we also incor- 
porate and reana lyze archival XMM- A^ewton and CXO 
observations (see IGavriil et al.l [2006t . for details on how 
these data were processed). Separate post-outburst op- 
tical and near-IR observations of IE 1048.1— 5937 are 
reported in detail elsewhere (Wa ng et al.ll2007l ). 

2.1. RXTE 

We have observ ed IE 1048.1—5937 reg ularly since 
1997 with RXTE ( Kasvi et al.|[200l IGK04I ). Our data 
were obtained using the Proportional Counter Array 
(PCA) on board RXTE which consists of five identical 
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and independent Xenon/Methane Proportional Counter 
Units (PCUs). We use our RXTE observations of 
IE 1048.1—5937 to look for the presence of bursts (see 
iGavriil et al.ll2004l for details), to look for pulse profile 
changes, to monitor its pulsed flux, and to monitor its 
frequency evolution using phase-coherent timing when 
possible. 

For the timing analysis, we created barycentered 
lightcurves in the 2—5.5 keV band with 31.25 ms time res- 
olution. As past monitoring has shown that it can be dif- 
ficult to maintain pulse phase coherence over timescales 
longer than a few weeks, observations of IE 1048.1—5937 
since 2002 are done three times per week, with the three 
observations carefully spaced so as to allow a phase- 
coherent analysis and a precise frequency measurement. 
Thus, for each observation, we fold at the pulse period 
determined via periodogram, cross-correlate the folded 
profiles with a high signal-to-noise template, and fit the 
resulting phases with a linear function whose slope pro- 
vides the average frequency. Frequencies determined in 
this way are shown in Figure[TJ As can be seen in this fig- 
ure, the 2004-2006 frequencies were actually much more 
stable than in the past; we therefore attempted a fully 
phase-coherent analysis as well. This will be described 
elsewhere (Dib et al., in preparation). Frequency deriva- 
tives, which are also displayed in Figure [1] and are from 
IGK04I and Dib et al., (in preparation), are clearly much 
more stable in the 2004-2006 interval than previously. 

This period of rotational stability was accompanied 
by X-ray pulsed flux stability and relative quiescence: 
the pulsed flux time series of IE 1048.1—5937 in the 
2—10 keV band is presented in Figure [TJ The pulsed 
flux was c alculated using a rn ethod similar to that de- 
scribed in IWoods et al.l (|2004l ) and is based on the rms 
of the folded profile, but without variance subtraction. 

The long period of rotational stability end ed in 2007 
March , when a large glitch was observed (jPib et al.l 
I2007H) . The glitch also signaled the end of the period 
of pulsed flux quiescence. The pulsed flux suddenly in- 
creased by a factor of ^^3 in the energy range 2—10 keV. 
The upper limit on the rise time of the pulsed fiux for 
this event is approximately one week. The peak pulsed 
flux reached by the source was ~10% larger than the 
peak flux reached during the largest of the two previ- 
ously observed flares, and the rise time was at least 4 
times smaller. On 2007 May 17, the date of the last 
X-ray imaging observation included in this paper, the 
pulsed flux had decreased by ^^10%. It is as yet difficult 
to compare the decay timescale of the pulsed flux of this 
new event to that of the previous flares however it will 
be possible in the near future. 

2.2. CXO and Swift 

X-ray imaging observations were carried out with the 
CXO and Swift telescopes. The date, total exposure 
time and resulting count rate for each region-filtered 
background-subtracted observation are listed in Table [T] 
and a detailed description of our spectral fitting is given 
in§[mi 

IE 1048.1-5937 was observed five times*^ in 2006 with 

* Although four simultaneous observations were originally 
planned, HST experienced technical difficulties during one, result- 
ing in a fifth CXO observation. 
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Fig. 2.— The evolution of IE 1048.1-5937's X-ray spectral prop- 
erties. For consistency, all the spectra from CXO, Swift, and XMM 
were fit jointly with a two-component absorbed blackbody plus 
power-law model, that produced a best-fit Nu = (0.97 ± 0.01) X 



errors only. Vertical dashed lines indicate the beginnings of the 
flux flares (see Fig. [T] caption) , while horizontal dotted lines indi- 
cate quiescent values, (a) Blackbody temperature kT. (6) Photon 
index F. (c) Ratio of the 2—10 keV fiux contribution from the 
blackbody and power law components, (d) Total unabsorbcd flux 
in 2—10 keV. (e) 2—10 keV rms pulsed fraction calculated using 
the method described in Woods et al. (2004). (/) For reference, we 
show the 2—10 keV rms pulsed fiux as observed by RXTE. 

CXO in quasi-equispaced intervals, and three times in 
2007 April. The observations were made with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) camera us- 
ing the S3 chip in continuous clocking (CC) mode. This 
mode generates a 1x1024 pixel image, which is read out 
every 2.8 ms; pile-up in this mode is not an issue. The 
subsequent analysis was based on the "Level 2" events 
files for which the event times are photon arrival times 
(as opposed to readout times for the "Level 1" events). 
Following the standard threads we then corrected the 
"Level 2" events for certain caveats related to ACIS CC- 
mode data^. We extracted a rectangular region centered 
on the source with a width typically of ~30 pixels. We 
estimated the background using rectangular regions on 
either side of the source which extended from ^10 pixels 
beyond the edge of the source region to ^10 pixels from 
the edge of the image. Following the standard threads, 
the source and background spectra were then extracted 
using CIAO^'^ v3.4. We grouped the source spectrum such 
that there were no fewer than 20 counts per bin after 
background subtraction. The response matrix function 
(rmf ) and area response function (arf) for each observa- 
tion were also generated using CIAQ and CALDB. We used 
Timed Exposure (TE) mode response matrices because 

^ http: //cxc .harvard. edu/ciao3 .4/why/ccmode .html 
http://cxc.harvard.edu/ciao 
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of the absence of available spectral calibration for CC 
mode. 

IE 1048.1-5937 was observed three times" in 2007 
with Swift using the X-ray Telescope (XRT). XRT 
has two observing modes: photon counting (PC) and 
windowed-timing (WT). PC mode provides a ^600x600 
pixel^ image with low (^2.5 s) time resolution, whereas 
WT mode provides a one-dimensional ~200-pixel wide 
image with high (^^0.74 ms) time resolution. XRT au- 
tomatically switches from WT to PC-mode when the 
count rate exceeds ~2 counts s^^. In our analysis we 
ignored the WT mode data as they were too short to 
be of use. We used the PC-mode event files from the 
standard pipeline; these were cleaned and barycentered. 
Using FTDDL xrtcentroid, we determined the centroid 
of the image, and input the cleaned, barycentered events 
into the command line interface xselect^^. We selected 
a circular source region with a 25 pixel radius centered on 
the centroid of the image. Since the source was slightly 
piled up, we excluded a 4 pixel radius circular region 
centered on the centroid of the image. We used 4 pixels 
because that is the recommended exclusion radius for a 
source with IE 1048.1— 5937's spectrum (L. Angelini, pri- 
vate communication). For the background, we selected 
events from an annulus with an inner radius of 50 pixels 
and outer radius of 130 pixels. We further filtered our 
events by selecting only those that had grades to 12. 
With xselect we then created source and background 
PHA files. The FTOOL xrtmkarf was used to create the 
arfs, and we input the source PHA file with the excluded 
center so that it would be corrected for pile-up. We used 
the file provided by CALDB appropriate for PC mode data 
events of grades to 12 to make the rmf. 

2.2.1. Global Spectral Fit 

We extracted spectra of IE 1048.1-5937 from all 
archival, monitoring, and ToO observations collected 
from CXO, XMM and Swift. A description of the re- 
cent data analysis is in ? I2.2|, and the a rchival anal- 
ysis was described in i Gavriil et al.l (|2006f) . Using the 
fitting package XSPEC^^ vl2.3.1 we modeled the X- 
ray spectra with a photoelectrically absorbed black- 
body plus power law. We used the phabs photoelec- 
tric absorption XSPEC model, which assumes th e so- 
lar abundance table of lAnders fc Grevessi (|1989( ) and 
uses the bcmc photoionization cross-sect ion table from 
iBalucinska-Church fc M cCammon (1992') with the new 
He cro ss-section from iBalucinska- Church fc McCammonI 
(|1998l ). We fit this model to all the observations si- 
multaneously, allowing the column density Nu to vary 
but with the sole constraint that it be the same for all 
observations. We restricted our fit to the 0.7-5 keV 
band. The total global fit had xl = 1-07 for v = 
4200 degrees of freedom (dof). Individually, the ob- 
servations were equally well modeled (except for a rel- 
atively high value in one observation, see § I2.2.2p . 
Our best-fit column density obtained from the global 
fit is Nh = (0.97 ± 0.01) X 10^^ c m^^, in agreement 
with lDurant fc van Kerkwijkl (j2006bl ). The results of our 

Two additional Swift observations taken contemporaneously 
were omitted from our analysis because their short exposures 
yielded prohibitively largo uncertainties. 

http://heasarc.gsfc.nasa.gov/docs/ sof tware/lheasof t/f too! 

http://xspec.gsfc.nasa.gov 



spectral fitting are listed in Table [T] and plotted in Fig- 
ure H 

2.2.2. Possible Spectral Feature 

All observations had spectra that were well modeled by 
a photoelectrically absorbed blackbody plus power law. 
However, the CXO observation on 2007 April 6 (obser- 
vation ID 7647), the first after the 2007 March event, 
had a relatively high — 1-33 {v — 353 dof). There 
is possible evidence for an absorption line at ~2.7 keV 
(see Fig. [3]). Adding a Gaussian fine improved the fit 
Ixl = 1-20 for V = 351 dof) with Ax^ = 48.3. With the 
line width fixed at 0.1 keV, we measure a line energy of 
2.73 ± 0.03 keV. Allowing the line width cr to be a free 
parameter, we obtain a — 0.16 ± 0.03 keV, a line energy 
of 2.74 ± 0.03 keV and xl = 1-19 for 351 dof. We find 
no evidence of a phase dependence for this possible line. 
A detailed phase-resolved spectroscopic analysis will be 
presented in a forthcoming paper. 

To test the significance of adding such a line, we per- 
formed the following simulation. We generated 10000 
simulated spectra by adding Poisson noise to a black- 
body plus power-law model spectrum having the same 
parameters as our best-fit model. We then determined 
the maximum change in x^ after adding a Gaussian line. 
To avoid local minima and ensure that we found the true 
minimum x^ for each simulation iteration, instead of fit- 
ting for the peak energy, we stepped through different 
line energies between 0.6 and 7.0 keV with a step size 
of 0.025 keV. The width of each line was held fixed at 
0.1 keV and its normalization was allowed to vary. In 
all the simulation iterations, none had a change in x^ 
greater than 48.3. The probability that this feature is due 
to random chance is <0.13%, accounting for the number 
of trials. We note however that with the addition of the 
line, the overall fit is still unacceptable; we speculate that 
other features may be present as well, though at marginal 
significance. For example, there is a small but intriguing 
feature at ~1.3 keV, which is approximately half the en- 
ergy of the line discussed above. The evidence for a line 
must be considered tentative because CC mode, which 
uses the TE mode response matrices, is not spectrally 
calibrated, and there exist calibration lines between 1.7 
and 3 keV in the TE response. However, this is true for 
the other CXO observations, including those in which we 
find no such features. 

2.2.3. Pulse Morphology and Pulsed Fraction Study 

Using our CXO and Swift observations we were also 
able to study the source's pulse morphology, flux and 
pulsed fraction. 

For the CXO observations, we barycentered our cor- 
rected Level 2 event lists and then extracted background- 
subtracted light curves in the different bands using the 
same source and background regions as for the spectral 
analysis. We folded these lightcurves at the optimal fre- 
quency as determined by a periodogram. The frequencies 
we obtained all agreed with the ephemerides determined 
with our contemporaneous RXTE data. Figure [4] (Ze/i) 
displays the CXO pulse profiles in the 1—3 and 3—10 keV 
bands. Notice that after the flare (last 3 panels in Fig- 
ure |4] left), the pulse profile changed from single to at 
s/Ssi§.fectff;BkiBi§}ffidi-i To study the pulse morphology evo- 
lution quantitatively we decomposed the pulse profiles 
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TABLE 1 

CXO, XMM AND Swift OBSERVING PARAMETERS AND RESULTS 



Date 


MJD 


Obs. ID 


Exposure 


Count rate kT'^ 


Unabs. Flux'^-t' 




Pulsed 








(s) 


(counts s^-"-) (keV) 


(lO^-"-^ erg cm^^ s^-"-) 




Fraction'* 


Archival XMM Observations 


2000 Dec 28 


51906 


0112780401 


4447 


1.461(18) 3.31(22) 0.656(17) 


4.8(4) 


0.6(2) 


0.76(2) 


2003 Jul 16 


52836 


0147860101 


41145 


3.429(9) 2.93(4) 0.599(6) 


11.5(1) 


1.21(5) 


0.44(4) 


2004 Jul 9 


53195 


0164570301 


32679 


0.985(6) 2.75(5) 0.542(15) 


9.6(1) 


2.1(2) 


0.561(7) 


Archival CXO Observations 


2004 Jul 10 


53196 


4653 


28860 


1.363(7) 2.87(7) 0.56(1) 


9.3(1) 


1.2(1) 


0.56(1) 


2004 Jul 15 


53201 


4654 


28085 


1.317(7) 2.87(7) 0.55(1) 


9.4(2) 


1.2(1) 


0.55(1) 


Monitoring CXO Observations 


2006 Feb 26 


53792 


6733 


22085 


1.072(8) 2.79(6) 0.51(1) 


7.1(1) 


2.0(2) 


0.63(1) 


2006 Apr 20 


53845 


6734 


20553 


1.157(9) 2.75(5) 0.50(1) 


7.1(2) 


2.0(2) 


0.61(1) 


2006 Jun 14 


53900 


6735 


22085 


0.985(9) 2.95(9) 0.55(1) 


6.3(1) 


1.4(2) 


0.64(1) 


2006 Jul 30 


53946 


7347 


22176 


0.821(8) 2.72(8) 0.51(1) 


6.8(2) 


1.2(2) 


0.70(1) 


2006 Sep 23 


54001 


6736 


22022 


1.006(7) 3.14(11) 0.56(1) 


5.1(1) 


1.0(2) 


0.64(1) 


ToO CXO Observations 


2007 Apr 6 


54196 


7647 


20079 


4.241(15) 2.30(4) 0.56(1) 


36.9(5) 


1.6(1) 


0.299(5) 


2007 Apr 16 


54206 


7648 


18387 


3.959(15) 2.27(4) 0.55(1) 


36.0(5) 


1.7(1) 


0.286(6) 


2007 Apr 28 


54218 


7649 


19082 


4.332(15) 2.12(4) 0.55(1) 


36.0(5) 


1.8(1) 


0.287(5) 


ToO Swift Observations 


2007 Apr 3 


54193 


30912001 


4905 


0.58(1) 2.7(2) 0.58(6) 


35(2) 


2.0(6) 


0.22(3) 


2007 Apr 13 


54203 


30912003 


4829 


0.55(1) 2.8(4) 0.62(4) 


32(2) 


0.7(3) 


0.17(3) 


2007 May 17 


54237 


30912009 


3626 


0.51(1) 2.3(4) 0.7(1) 


35(3) 


2(1) 


0.22(3) 



Note. — All quoted err o rs represent la unce rtain ties. The arch i val X MM and CXO observations have been previously reported 
bv lMereghetti et al.l l l200i) . iTiengo et al.l (I^OOST) . and lGavriil et al.l ||2006| ). 

^ These results are from a simultaneous spectral fit to all the data using as a model a photoelectrically absorbed blackbody with 
temperature kT and a power law with photon index F. In the simultaneous fit the column density, Njj, was constrained to be the 
same for all observations. The resultant Njj from the global fit was Njj = (0.97 ± 0.01) X 10^^ cm~^; see the text for details.'' Total 
unabsorbed fiux in the 2—10 keV band.'^ Ratio of the 2—10 keV flux in the power-law component, -FpL, to the 2—10 keV flux in the 
blackbody component, Fbb.'^ 2—10 keV pulsed fraction, using the definition of pulsed fraction described in lWoods et al.l I I2004I) . 
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Fig. 3.— IE 1048.1-5937's X-ray spectrum on 2007 April 6, 
as observed by CXO. Top: The best-fit photoelectrically absorbed 
blackbody plus power law is plotted as a solid line; see text for 
details on the analysis. The quality of the fit was low, with = 
1.33 for V = 353 dof. Bottom: Residuals after subtracting the 
best-fit model. Notice a possible spectral feature at ~2.7 keV. 

into their Fourier components. In Figure 3] (right) we 
plot the Fourier components in terms of their ratio to 
the fundamental (the n — Q and fundamental are ex- 



cluded from the plot). Notice that before the flare the 
profile consisted only of the fundamental and an n = 2 
Fourier component. After the flare, the profile exhibited 
significant Fourier components of order n = 2, 3, 5. Inter- 
estingly, despite the drastic changes in the other Fourier 
components, the n — A component remains consistent 
with zero throughout the flare. 

Similarly, for the Swift observations we used the 
barycentered, cleaned, and pointed PC-mode data 
events. The XRT WT data consisted of too few pulses 
to be useful. Swift PC-mode data have only 2.5-s reso- 
lution, providing fewer than 3 pulse phase bins for this 
6.45-s pulsar. We attempted to improve on the coarse 
sampling by utilizing advanced "bin splitting" techniques 
(G. Israel, private communication), however we find only 
marginal evi dence for a pulse profile change from these 
data (but see lCampana fc Israelll2007b . 

Using our folded pulse profiles we measured root-mean- 
square (rm s) pulsed fra c tions^ "^ using the method de- 
scribed in IWoods et all (|2004D . Our 2-10 keV rms 
pulsed fractions are listed in Table [1] and plotted in Fig- 
ures [J and The XMM pulsed fractions, taken from 
iGavriil et alT(l200l) . were extracted in the same fashion. 
Because of the low time resolution of the Swift data, the 
rms pulsed fraction would be artificially reduced if the 
pulse contained as much harmonic content as the CXO 
data during the flare. Thus, for each Swift observation we 
simulated, given its observed pulsed fraction and mean 

For a simple signal of the form c((/i) = Asin(27r(?!>)-|-iJ, the rms 

^ ^ , and the peak-to-peak 



V2 B 



pulsed fraction is given by PFrms = 

pulsed fraction is given by PFpp = We opted to use the rms 
with Fourier filtering (Archibald et al., in preparation) because it is 
a more robust estimator of the pulsed fraction than peak-to-peak. 
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Fig. 4. — Left: 1—3 keV {blue) and 3—10 keV {red) normalized pulsed profiles as observed by CXO. The profiles are phase aligned, have 
had their minimum bins subtracted, and are normalized such that the area underneath them is unity. Note that the 2007 event occured 
around MJD 54183. Right: 1—3 keV {blue) and 3—10 keV {red) Fourier decomposition of the pulse profiles. The powers P„ are plotted 
in terms of their ratio to the fundamental Pi. The n = and the n = 1 (fundamental) components are excluded. Notice how the pulse 
profile has additional structure after the glitch/fiux enhancement. 

count rate, what its true pulsed fraction would be as- 
suming a pulse profile with as much harmonic content as 
the nearest CXO observation. On average we found a 
reduction of no more than ~73%. Note in Figure [5] the 
observed tight c orrelation between t otal flux and pulsed 
fraction (see also lTiengo et al.1l2005() which we discuss in 
detail in § [5?^ 

2.3. HST 

We observed IE 1048.1-5937 using HST simuhane- 
ously with CXO on four occasions in 2006: February 26, 
April 20, July 30, and September 23 (Program ID 10761). 
Observations were made with the Near Infrared Camera 
and Multi-Object Spectrometer (NICMOS) instrument, 
a 256 X 256 square pixel HgCdTe detector. Combined 
with camera 3, it provided a focus ratio of //17, a field 
of view of 51' X 51', and a plate scale of 0'.'2 pixel^^. The 
detector was read out in MULTIACCUM mode, and a 
spiral dither pattern with 5" spacings was applied. We 
observed using the filter FllOW (similar to the ground 
based J-band filter) on two occasions and F160W (sim- 
ilar to i?-band) on four occasions. The average FWHM 
of the point spread function (PSF) was ~ 0'.'4 for all 
NICMOS observations. Observing parameters are sum- 
marized in Table [21 

The data underwent "On-The-Fly-Reprocessing" in- 
volving standard pipeline routines before being retrieved 
from the HST Archive. The calibration routine calnica 
performed basic data reduction steps, and calnicb pro- 
duced a single final image for each observation, com- 
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prior to that are from previous literature llWang 
\2Q0% HsraeTe t al. 2002; Durant ' " 



53000 54000 
Modified Julian Days 

Fig. 5.— The evolution of IE 1048.1-5937's near-IR brightness. 
VLT and HST results from this work fall after 200 4, while data 

„ Chakrabartvl 

van Kerk wiik 20051. Wc also 
show the corresponding magnitude on the right axis; however the 
translation is not exact due to HST's unique filters. Vertical 
dashed lines indicate the beginnings of the flux flares (see Fig. [Tl 
caption), while horizontal dotted lines indicate quiescent values, 
(a) "J-band" flux i/F,^. {b) "iJ-band" flux uF,y. The detection 
made from combining three HST observations (see S I2.3| l is shown 
in grey, (c) "A",<;-b and" flux vFu. Res ults from 2007 are described 
in detail elsewhere l lWane et al.ll2007l ). {d) For reference, we show 
the 2—10 keV rms pulsed flux as observed by RXTE. 
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bining the dit hered mosaic input images. DAOPHOT 
(| Stetson! [T987h for IRAF v2.12.2 was used to perform 
PSF photometry. Absolute cahbration of our measure- 
ments was done using the NICMOS Photometric Key- 
words'^ for which the uncertainties are beheved to be 
less than 5%; we found the aperture correction to corre- 
spond with the calibration keywords by simulating PSFs 
with the Tiny Tim'^ software package. 

In order to determine limiting instrumental magni- 
tudes in each of the observations, we performed PSF 
photometry on a set of a hundred trial images to which 
a single artificial star was added. This artificial star was 
placed in a (X'S box in a blank region near the nominal 
position of the IR counterpart. By varying the bright- 
ness of the star, we define the 3cr detection limit at the 
instrumental magnitude for which the PSF photometry 
reco vers the artificial star with an uncertainy of 0.3 mag 
(see iHulleman et~aIll2000D . 

At the position of I E 104 8.1—5937, as measured by 
IWang fc Chakrabartvl (|2002t l in Ks-hand with the Mag- 
ellan telescope, we find one point source in the F160W 
image from 2006 February 26 only. To confirm the po- 
sitional coincidence of this source with the Magellan ob- 
ject, we astrometrically tied = 31 nearby field stars 
in our HST image to those from the Magellan image. 
A - O'.'OS offset was found for the IE 1048.1-5937 
candidate. The Magellan source's relative positional 
uncertainty, given its 0'.'4 FWHM radius, is roughly 
FWHM/V^ - 0'.'07. Therefore, we consider this HST 
object the counterpart of IE 1048.1-5937. 

Thus, we measure our lone detection of 
IE 1048.1-5937 to have mpieow = 22.70 ± 0.14 mag; 
uncertainties indicate the nominal errors determined by 
DAOPHOT and incorporate PSF fitting uncertainties. 
Limiting magnitudes for the three latter HST F160W 
observations in which no counterpart was seen are pre- 
sented in Table [21 To deepen the sensitivity, geometric 
transforms were computed to register and combine the 
three non-detection datasets. In the combined data, 
we tentatively detect an extremely faint counterpart 
of brightness mpi^ow = 23.42 ± 0.20 mag, where the 
3ct upper hmit is > 23.5 mag. We find no FllOW 
counterpart down to a 3a upper limit of mpiiow 
> 25.0 mag and > 24.8 mag, on February 26 and April 
20, respectively. Results are plotted in Figure O 

To examine how this compares with past results, we 
also calibra ted our measurements u s ing t he H magni- 
tudes from iDurant fc van KerkwijS ()2005l ). A relative 
offset is found between HST's F160W filter and ground 
based i/-band, as a function oi J — H color. We know 
the color of IE 1048.1—5937 onl y during a brightened 
state {J -H K 0.9 in April 2002; IWang fc Chakrabartvl 
|2002| ). and suspect that it is not constant; therefore, this 
is not a perfect comparison, but provides a rough esti- 
mate. Based on the magnitudes of seven nearby field 
objects, we measure H m 22.73 mag on February 26, 
which is very near what we find for mpimw- 

2.4. VLT 

We have analyzed archival observations of 
IE 1048.1-5937 which were obtained with NAOS- 

http: //www. stsci . edu/hst/nicmos/perf ormance/photometry 
http: //www. stsci . edu/sof tware/tinytim/tinytim.html 



CONIC A (NACO), the near-IR adaptive optics (AO) 
instrument at VLT of the European Southern Observa- 
tory. The source was observed for 900 s in i^Tg-band using 
the S27 camera (28" x 28" field-of-view and 27 maspix"' 
pixel scale) on four occasions (A pril 11, April 29, May 
27 an d June 6) in 2005. Star A (jWang fc Chakrabartvl 
I2002D . located 9'.'5 from the AXP, was used as the 
wavefront-sensing (WFS) star. Due to the relative 
faintness of the WFS star [r' = 16.5 mag), the AO 
correction was only partial. 

The observations, either dithered sets of 5 2 x 90 s inte- 
grations (on the first night) or of 10 1 x 90s integrations, 
were corrected for dark current, fiat-fielded (using aver- 
ages derived from the science frames) and corrected for 
variations in the sky (using median averages of the fiat- 
fielded science frames). Finally, the reduced science im- 
ages were aligned using integer pixel offsets and median- 
combined to create an average image for each night. 

Using DAOPHOT, we tried PSF photometry on each 
of the four averaged images, using different models for 
the PSF. Because the AO correction is only partial, we 
found that the PSF is reasonably well modeled by a Mof- 
fat function of exponent 2.5 in combination with a look- 
up table. The width of the PSF that was determined for 
each observation was 0"34 (April 11), 0"43 (April 29), 
0"24 (May 27) and 0"20 (June 6). Only in the observa- 
tion of May 27 did we detect the near-IR counterpart to 
IE 1048.1-5937. 

We calibrated the instrumental magnitudes and mag- 
nitude limits using stars X2 to X8 and A and B against 
the calibrated Ks magnitudes of Wang & Chakrabarty 
(2002). The rms uncertainty in the calibration is about 
0.2 mag. The near-IR counterpart to IE 1048.1-5937 
was detected at = 21.0 ± 0.3 on 2005 May 27; the 
3(7 upper limits for April 11, April 29, May 27 and June 
6 are > 20.3, > 20.6, > 21.2 and > 21.4, 
respectively. Limiting magnitudes were determined with 
the same method described in § 12.31 Results are shown 
in Figure [5] 

Note that post-2007 event optical and near-IR obser- 
vations of IE 1048.1-5937 were obtained with VLT and 
Magellan (IWang et al.ll2007t l!srael et al.ll2007D . Prelim- 
inary results of the Ks-ha,nd observations are shown in 
Figure m and details on som e of the observations can be 
found in lWang et al! (|2007D . 

3. DISCUSSION 

Explaining the origin of the different types of AXP 
variability across the electromagnetic spectrum is a ma- 
jor challenge for the magnetar model. Many fundamental 
questions remain open. In the case of IE 1048.1—5937, 
how are timing instabilities related to radiative variabil- 
ity? AXP spin-up glitches sometimes are accompanied 
by major radiative events, as in the 2002 outburst of 
AXP IE 2259-1-586 ()Kaspi et al.ll2003l ). and som etimes 
not, a s in the large glitch of IE 1841-045 (Di b et al] 
l2007al) . Moreover, are the gradual large-scale changes 
seen in the 2001-2002 fiares of IE 1048.1-5937 of a sim- 
ilar origin to AXP o utbursts like that o f IE 2259-1-586 
or XTE J1810-197 (jlbrahim et al.ll2004D . and how can 
they be accounted for by magnetar models? Also, what 
is the nature of the optical/IR emission observed in five 
out of eight confirmed AXPs? Suggested origins include 
coherent or curvature emission in the plasma magne- 
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TABLE 2 

HST AND VLT OBSERVING PARAMETERS AND RESULTS 



Date 


MJD 


Exposure 


Filter^' 


Limiting 


Detection 








(s) 




Magnitude 




(lO"-""^ ergs s"-"- cm^'^) 


VLT Observations 


2005 Apr 11 


52324 


900 


Ks 


>20.3 




<7.0 


2005 Apr 29 


53471 


900 


Ks 


>20.6 




<5.3 


2005 May 27 


53489 


900 


Ks 


>21.2 


21.0(3) 


<3.1; 3.7±1.0 


2005 Jun 6 


53527 


900 


Ks 


>21.4 




<2.6 


HST Observations 


2006 Feb 26 


53792 


2637 


FllOW 


>25.0 




<0.37 


2006 Feb 26 


53792 


1196 


F160W 


>23.7 


22.70(14) 


<0.49; 1.3±0.2 


2006 Apr 20 


53845 


2637 


FllOW 


>24.8 




<0.45 


2006 Apr 20 


53845 


1196 


F160W 


>23.1 




<0.89 


2006 Jul 30 


53946 


1037 


F160W 


>23.3 




<0.74 


2006 Sep 23 


54001 


1037 


F160W 


>22.9 




<1.0 



Note. — Magnitude uncertainties reflect errors determined by DAOPHOT; upper limits are 
3(7 limits, as defined in the text. Measured quantities are "observed" , ie. have not been corrected 
for reddening/extinction effects. 

Although the filter system of HST does not precisely match standard ground-based near-IR 
filters, note that FllOW » J-band and F160W » H-ha.nd}' To convert standard IR magnitudes 
to fl ux, we take as the Ks = mag zero point vF^, = 9.28 X lO"'^ ergs s~^ cm~^, derived from 
I I2000I) . HST fluxes are determined from NICMOS Photometric Keywords. 



tosph ere (|Eichler et al.ll200^ iBeloborodov fc ThompsorJ 
|2007| ). but so far this has not be confirmed. Are the 
low- and high-energy emission mechanisms intimately 
connected (Hcyl & Hernquist 2005)? Correlated post- 
burst flux decay in the X-ray and n ear-IR regimes h as 
been seen in at least one other AXP ()Tam et al.ll2004D . 

Given the many observational properties that have 
been characterized for IE 1048.1—5937 in this work and 
for other AXPs in similar studies, which property will 
emerge as being the most constraining of physics is hard 
to know. The most promising behaviors are those which 
show clear correlations with others, or those that are 
common to many or all AXPs and SGRs. In this section, 
we discuss the behavior of IE 1048.1—5937, and focus on 
the phenomena that are potentially the most useful for 
testing the magnetar or other competing models. We 
first consider the quiescent phase we have observed in 
2004-2006, and its implications for AXP models, then 
subsequently discuss the source's return to activity in 
2007 March. 

3.1. The 2004-2006 Quiescent Phase 
3.1.1. X-ray Flux and Spectrum 

The X-ray flaring observed pre-2004 contrasts strongly 
with the stable pulsed fluxes we observed in 2004-2006 
(Fig. [T]c). Clearly the source's stable state is also its 
faintest. The X-ray spectrum during this quiescent phase 
is also fairly constant, with kT evolving post-flare, on a 
time scale of several years. Interestingly, kT does not 
return to its pre-flare value as measured in 2000 Decem- 
ber (see Fig. [2]3). Similarly, the photon index (Fig. [2j)) 
decreased in 2004-2006, away from its softer pre-flare 
value. Meanwhile the pulsed fraction (Figs. [T]e and [Us) 
was slowly rising in 2004-2006, as the source grew fainter, 
as if, unlike kT and P, it was slowly recovering to the pre- 
flare value. Thus, overall, the quiescent period is char- 
acterized by slow evolution, on timescales of years, in 
which the source flux and pulsed fraction slowly relaxed 
back to their pre-flare values, while the source spectrum 
varied significantly (Fig. [Uc) and did not relax back to 
its pre-flare state. This is suggestive of a fixed energy 



loss rate in quiescence, though perhaps with a different 
magnetospheric current configuration, which impacts the 
surface emission via return currents ([Thompson et aH 
,2002) . We note, on the other hand, that the 2006 and 
pre-2001 pulse pro file were ver y similar if not identica l 
(|Kaspi et al]l2001h . and that the [Thompson et al.l ([2002[ ) 
simple prediction for the decay time of the sort of magne- 
tospheric twist required for IE 1048.1—5937 is an order 
of magnitude too large for standard parameters. The 
evolution could also be purely due to the thermal com- 
ponent ([Ozel fc Guveij[2007f) . 

3.1.2. Timing Stability 

Simultaneous with the slow evolution of the X-ray flux, 
pulsed fraction, and spectrum, the pulsar's rotational be- 
havior clearly stabilized, with the source spinning down 
relatively smoothly at a value close to the long-term aver- 
age in 2004-2006 (Fig.[Th). Whereas previously this AXP 
distinguished itself from others by defying attempts at 
phase-coherent timing and exhibiting large torque vari- 
ations (Fig. [Ij)), in 2004-2006, IE 1048.1-5937 resem- 
bled, from a timing point of view, other AXPs which 
are r elatively st able rotators, at least when not glitching 
(e.g. [Kaspi eTal . 1999; Dib et al. 2007a). This behavior 
demonstrates a clear relationship between timing and ra- 
diative properties, and that the "noise" seen in the fre- 
quency derivative during 2001-2004 is likely physically 
different from the "timing noise" observed ubiquitously 
in radio pulsars and in otherwise radiatively stable AXPs. 

In the context of the magnetar model, magnetospheric 
activity can account for both torque and X-ray fLux vari- 
ability, although the former is most sensitive to currents 
anchored closest to the magnetic poles, so that only a 
broad correlation betw een flux and torque is expected 
([Thompson et al.ll2002D and in deed is seen when consid - 
ering many AXPs and SGRs ([Marsden fc WhitdlMilh . 
That the torque and luminosity in IE 1048.1-5937 
do not vary simultaneously or in a clearly correlated 
way is thus not necessarily problematic if the magne- 
tospheric current configuration is changing, although it 
does demonstrate that the magnetar model in this par- 
ticular regard is not strongly predictive. 
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Although the magnetar model is favored, it has also 
been suggested that accretion from a fossil debris disk 
could explain AXP spin characteristics (jChatteriee et al.l 
UPOO; Alpar 20d l|) as well as all asp ects of the broad band 
emission (Ertan & Caliskan"2006). While it IS now ev- 
ident that accretion alone cannot be responsible for all 
observed properties, most notably the energetic X-ray 
bursts, a "hybrid" model has been invoked that puts 
a thin debris disk around a highly magnetized pulsar 
(|Eksi fc Alpaill2003 f). In this hybrid case, both the per- 
sistent luminosity and pulsar spin down are related to 
the mass transfer rate, M. According to IGK04I . who 
compare IE 1048.1— 5937's torque changes with changes 
in RXTE pulsed flux, the scale of their variability does 
not obey an expected relationship for an accreting pul- 
sar undergoing spin-down, thus presenting a challenge 
to models of fossil disk accretion. It could be argued, 
however, that since only the pulsed flux Px was being 
monitored, the total X-ray luminosity Lx is still an un- 
known quantity. We address this by using the quanti- 
tative correlation we establish between pulsed fraction 
Vj^ and total unabsorbed flux Fx (discussed in § 13.2.21 
and shown in Fig. [S]) to simulate a well sampled set of 
phase-averaged flux data as a function of RXTE pulsed 
flux Px , given that Vjr is simply Px / Fx ■ The resulting 
"new" total flux. Fx, is shown in Figure [TJi, along with 
v in Figure [T]!) for comparison. At its most variable in 
2002-2004, the absolute value of v changed by a factor 
of >10 in less than a year, while the maximal change in 
total unabsorbed flux is a factor of ^6 from peak to qui- 
escence. Note that these variations are not simultaneous. 
For a pulsar experiencing a spin-down torque due to mass 
accretion while in a quasi-equilibrium "tracking" phase 
(i.e. the AXP phase), a strong correlation described by 
Lx oc \v\''/^ can be derived from lChatterjee et al.l ()200d . 
equation 3), where the magnitudes of torque and v are 
proportional, the radius of the magnetosphere is defined 
by the Alfen radius which is also a function of M, and 
Lx oc M. Thus, a factor of >10 change in v should 
be reflected by a factor of >200 simultaneous change in 
Lx, and thus Fx, if this model is correct. Clearly this is 
not observed, rendering this particular accretion scenario 
unlikely. This point is further emphasized by the signifi- 
cant time offset between the changes in torque and flux. 
OveraU, the fact that IE 1048.1-5937's long history of 
highly irregular spin-down is mirrored by trends in its 
X-ray emission only in a broad, rather than strict, sense 
strongly suggests that active accretion is not happening 
in this case. 

3.1.3. Near-IR Quiescence 

Our original purpose for making simultaneous 
CXO / HST observations was to compare low-level near- 
IR and X-ray fiux and spectral changes, with each 
observed with the same instrument, in order to look 
for correlations. Given that the near-IR source had 
faded considerably in the HST observations, render- 
ing it only marginally detected in 2006, this was not 
possible. Compared to the handful of near-IR detec- 
tions made throughout 2002-2003 (see Fig. [5]) it is clear 
that 2005-2006 marked a period of near-IR quiescence in 
IE 1048.1-5937. At the time of the 2005 VLT observa- 
tions in which the AXP was faintly seen in iiTg-band, its 



flux was consistent with the last detection in mid-2003. 
With HST in 2006 we find that its H and J magnitudes 
have dropped lower than ever before observed. There is 
marginal evidence for near-IR flux variability during qui- 
escence on comparably short timescales to the variabil- 
ity in X-ray spectral parameters, although closely spaced 
near-IR monitoring observations are required to confirm 
this. 

That the near-IR faded roughly in concert with the 
X-ray fiux is notable, suggesting a correlation similar 
to that seen in AXP IE 2259+586 (jTam et all [2004 ). 
This is discussed further below. However such a cor- 
relation can be argued to be expected in both the 
magnetar and disk models. In the magnetar model it 
would imply that the near-IR emission is magnetospheric 



tospheric configuration varies ( 


Hevl & Hernouistl 120051: 


iBeloborodov & Thomosonl 12007 


). In the disk model, a 



correlation between near-IR emission and X-ray fiux is 
naturally expected since the putative disk is heated via 
X-ray illumination (Ertan et al. 2006; Wang e t al. 2006) , 
however why the torque should have varied strongly and 
non-simultaneously is a puzzle (see above). 

In any case, the concurrence of the near-IR and X-ray 
quiescence and the timing stability is unlikely to be an 
accident. Similarly, the presence of the two flares and 
the subsequent strong torque changes are also not likely 
to be by chance. 

3.2. The 2007 March Event and its Aftermath 

IE 1048.1-5937 reactivated in 2007 March, when a 
sudden pulsed flux increase was observed (Fig. [T]c) , ac- 
companied by a large spin-up glitch, the details of which 
will appear elsewhere (Dib et al., in preparation). Our 
follow-up X-ray observations with CXO and Swift show 
that the spectral parameters, pulse shape, total flux, and 
pulsed fraction are drastically changed since this flare oc- 
curred (Figs.[2]and|3]). In particular, the total X-ray flux 
was >7 times greater in the 2—10 keV band compared 
with quiescence, while the pulsed fraction decreased from 
~75% to ~20%. Compared to the previous long-term 
flares from IE 1048.1—5937, the onset of this event took 
place more quickly by a factor of >4, the peak pulsed flux 
is ^10% greater, and the "high" state appears to be last- 
ing longer. Optical and ne ar-IR follow-up w i th the Mag- 
ellan Telescope and VLT (jWang et al.ll2007t [Israel et all 
|2007| ) reveal an incr ease in both I an d i^s-band of ~1.3 
mag. Very recently, iRea et aH ()2007[ ) reported on XMM 
ToO observations, and revealed that nearly three months 
after the flare onset, the total X-ray flux has decreased 
slightly but is still ~5 times brighter than in quiescence. 

3.2.1. Near-IR Enhancement 

Previously, an anti-correlation between near-IR and 
X-ray flux in this AXP was suggested, based on the 
correspondence of the highest near-IR detection with a 
low p oint in X-ray pulsed flux^^ (jPurant fc van KerkwiikI 
|2005| ). Given the new measurements from both before 
and after the recent event, such an anti-correlation is 
highly questionable. The significant rise in near-IR flux 



Note that an error in Figure 4 of IDurant &: van Kerkwiifl 
l|2005l ) places the brightest Ks flux on the wrong date; we have 
corrected this in our Figure |5j;. 
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now appears correlated with the most recent X-ray flare. 
Interestingly, in 2002, a high near-IR detection was be- 
tween the two X-ray flares; this suggests that the near-IR 
emission is not exactly correlated with the X-rays, but 
rather an enhancement can precede, follow, or last longer 
than an X-ray enhancement. A near-IR brightening pre- 
ceding an X-ray enhancement would rule out illuminated 
disk models for this source; frequent mid- or near-IR ob- 
servations are required to check this. As discussed above, 
AXP IE 2259-1-586 also demonstrated strongly correlated 
near-IR and X-ray flux d ecay following its large 2002 
outburst (|Tam et al.ll2004( l. but in contrast, 4U 0142-^61 
was highly variable in the near-I R despite X-ray stability 
(jPurant fc van Kerkwiikll2006cD . and has shown no ev- 
idence of near-IR changes coinciding with X-ray bursts 
(Gonzalez et al., in preparation). Note that the claim 
of correlated flux decay in the case of XTE J1810— 19 7 
(|Rea et al.|[2004D is under dispute (|Camilo et al.ll2007D . 
Such inconsistent behaviour is puzzling for both the mag- 
netar and disk models. 

3.2.2. Correlation between Pulsed Fraction and Total Flux 

One conclusion to be drawn from Figures [Ij/ and [l]s 
is that IE 1048.1— 5937's pulsed fraction Vj^ is strongly 
anti-correlated wi th the total X -ray flux Fx , a tr end al- 
ready noticed by i Tiengo et al.l (2005) and Gavrii rit al.l 
([2006) . Now, with our larger sample of data covering 
a wider dynamic range, we can quantify this correla- 
tion: we find a power-law dependence of Vj^ oc 
where n = -0.46 ± 0.02 (see Fig. O, with Fx in the 
2—10 keV range. In fitting, we set the index as a free 
parameter, ie. not tied to any model; therefore, some 
scatter is expected, as evident in the large xl — 9.36 
produced by the fit {i/ — 14 dof). However, in find- 
ing the uncertainty on n, we have scaled the Vj^ data 
uncertainties to force = 1, effectively scaling up ct„ 
as well. A similar correlation has bee n proposed in an - 
other AXP, IRXS J170849. 0-400910 (iDib et al.ll20"07ah . 
although in that case it is such that the pulsed flux re- 
mains ne arly constant in the pr esence of >50% total flux 
changes (jCampana et al.|[2007[ ). 

How could an increase in phase-averaged flux be met 
with a simultaneous decrease in pulsed fraction? A grow- 
ing hot spot on the neutron-star surface, a result of ei- 
ther a changing magnitude or configuration of returning 
magnetospheric currents or internal processes, could at 
least in principle produce such an effect, if the initial hot 
spot size were large enough that its size as viewed from a 
distant observer were a significant fraction of the stellar 
surface. We note that lOzel &: Guveil (|2007[ ) suggest that 
magnetar afterglows may be dominated by surface ther- 
mal, rather than magnetospheric, changes; detailed simu- 
lations of AXP pulsations using gravitational light bend- 
ing and appropria te radiation beaming functions (e.g. 
iDedeo et a l."2001') in addition to magnetosp heric scat- 
tering ([Thompson et al.l 120021 : [Fernandez fc ThompsonI 
are needed to see if the observed correlation can 
be reproduced. 

3.2.3. Correlation between Hardness and Total Flux 

A pre diction made by the twisted magnetosphere 
model of Tho mpson et all (|2002f ) is that enhanced emis- 
sion should also be spectrally harder, since both are asso- 
ciated with larger magnetospheric twist angles. Indeed, 



we observed a correlation between the total X-ray flux 
and spectral hardness: see Figure [6l Similar behaviour 
has also been r e ported in IRXS J17 0849.0-400910 
()Rea et al.ll2005al: ICampana et al.ll2007D . This can be 
interpreted as a confirmation of an important magnetar 
model prediction, namely that both the spectral hard- 
ness and the total flux should increase for an increas - 
ing magnetospheric twist angle (jThompson et al.ll2002D . 
Hardening of the spectrum might also be an expected 
effect of increased luminosity if a large injection of ther- 
mal seed photons is repeatedly up-scattered due to reso- 
nant cyclotron scattering feedback processes in the mag- 
netosphere, thereby shifting the photon energies higher. 
Recently lOzel fc Guveil (|2007l ) have suggested that the 
hardness-intensity correlation in magnetars is a result of 
purely surface thermal changes, with no change in the 
magnetospheric configuration. It would be interesting to 
apply their model to the IE 1048.1—5937 data as they 
span a much greater dynamic range in flux than did the 
sources studied bv lOzel fc Guveil (|2007t) . This is beyond 
the scope of our paper, however. 

3.2.4. 2007 Pulse Profile Changes 

As seen in the CXO data obtained after 2007 March, 
the pulse profile of IE 1048.1—5937 changed abruptly 
after the X-ray flux enhancement and glitch. Several 
new harmonics are clearly visible, and are present in 
a largely energy- independent way (Fig. [4|). This extra 
power in the higher harmonics clearly signals a com- 
plication of the surface and/or magnetospheric config- 
uration. This change contrasts with that seen after the 
2002 outburst of IE 2259-1-586, which involved mainly an 
exchange of powers between the fundamental and first 
harmonic, with higher harmon ics remaining relatively 
unchanged () Woods et al.ll200^ . Also, the change con- 
trasts with that seen during the 1998 SGR 1900-1-14 giant 
flare, in which the SGR pulse profile simp lified greatly 
immediately post-flare (jGogiis et al.l I2002D . However, 
the opposite oc cured followed the 2004 giant flare from 
SGR 1806-20 (jPalmer et al.ll2005( l. when a sinusoidal- 
to-coniplex evolution in pulse morphology took place, 
similar to that of IE 1048.1-5937. This diversity of 
behaviors in magnetars post-outburst is problematic for 
understanding the underlying physics; indeed it indicates 
a wide variety of phenomenon phase space, which itself 
must be explained by models. The peculiar suppression 
of the 4th harmonic in the post-glitch IE 1048.1-5937 
profile is particularly puzzling to us; perhaps it indicates 
an important and unchangeable symmetry in the com- 
bined emission and viewing geometries, or perhaps it is 
a chance occurence and will not occur again in future 
events. 

3.2.5. 2007 April Spectral Feature? 

In the CXO observation made on 2007 April 6, the 
one closest to and just after the 2007 March event, 
the X-ray spectrum of IE 1048.1—5937 could not be 
well fit by a simple continuum model, in contrast with 
the other datasets. We identified a possible absorp- 
tion line near 2.7 keV, though we note that the spec- 
tral model including this line is still not a good fit to 
the data. This suggests that other low-level features 
may be present as well, including a possible feature at 



11 



o.t 



C 0.6 
o 

o 

!-, 

-a 

Oh 



0.2 




0.7 



0.6 



n3 



0.5 - 



cxo 

XMM 
Swift 



0.4 




10 
Unabs. 



2x10-11 3x10- 
2-10 keV Flux (erg/s/cm^) 



10-11 2x10-11 3x10-11 
Unabs. 2-10 keV Flux (erg/s/cm^) 



Fig. 6. — X-ray flux dependent properties of IE 1048.1—5937, from 6 years of observations with CXO (squares), XMM (triangles) , and 
Swift (circles). Left: Pulsed fraction vs. total unabsorbed flux, both in 2—10 keV. The solid black line shows the best-fit power law that 
describes the correlation between the two quantities. Right: Hardness ratio vs. total unabsorbed 2—10 keV flux. The hardness ratio is 
deflned as (S — H)/(S + H) where S and H are the phase-averaged count rates in the 1—3 and 3—10 keV bands, respectively. 



roughly half the above line's energy, ~1.3 keV. The pu- 
tative 2.7 keV line joins the growing list of puzzling 
and sometim es statistically marginal AXP/SGR spec- 
tral features (jlbrahim et al.ll2002l [200l iRea et al.ll2003l . 
l2005aD . Examples of some feature detections that are 
clearly statistically significant are, interest ingly, repeat- 
edly near -13 keV, e.g. in IE 10 48.1-5937 ("Gavrii l et all 
[2002. 2006.) . in XTE ,11810-197 tWoods et al. 20051) . and 
very recently AXP 4U 01424-61 (|Gavriil et ai.ll2007l . and 
Gavriil et al., in preparation). The feature we see in 
the 2007 April data for IE 1048.1-5937 are unhke any 
of these. It is perhaps more reminiscent of the absorp- 
tion features seen in some "isolated neutron st ars" (e.g. 
iBignami et "all 120031: Ivan Kerkwiik et al.ll2004D . and, if 
similarly interpreted as an electron cyclotron line, im- 
plies a magnetic field of — 3 x 10^^ G, or — 6 x 10^^ G 
for a proton cyclotron line. The latter value is of course 
more in line with the field expected near the surface of a 
magnetar. Ho wever, the feature could also be an atomic 
transition fe.g. iMori fc Holl2007f ): it is difficult to know 
given the present data. Regardless, the spectrum 10 days 
later, on 2007 April 16, showed no such features, indicat- 
ing that whatever their origin, it had subsided, and on a 
shorter timescale than that of any profile or ffux relax- 
ation. 

3.2.6. Comparison of the 2007 IE 1048.1-5937 Event with 
Other Similar Events 

The 2001-2002 flares observed in IE 1048.1-5937 are 
unique in th at th ey had very well resolved rise times of 
many weeks (GK04j). The 2007 April event, by contrast, 
involved an X-ray pulsed flux enhancement that rose at 
least 4 times faster. In this sense it is perha ps more sim- 
ilar t o the 2002 outburst of IE 2259-^586 (jKaspi et all 
|2003( ). although in that case a brief rise could have been 
missed by the sparse monitoring. During that outburst, 
IE 2259-f 586 exhibited an order-of-magnitude increase 
in pulsed and phase-averaged flux, as well as numerous 
other changes including (but not limited to) ^^80 SGR- 



hke bursts. Similarly to IE 1048.1-5937, both the to- 
tal and pulsed flux increased during the IE 2259+586 
outburst, and the pulsed fraction d ecreased, going from 
23% in quiescence down to 15% (|Woods et al.l l2004f) . 
However, the pulsed fraction was not as tightly cor- 
related to the pulse d flux: the puls ed fraction recov- 
ered within ^--^3 days (|Zhu et al.ll2007f) . whereas the flux 
took several m onths to return to its pre-outburst value 
(|Woods et al.l [2004) . In terms of energetics, the total 
energy in 2—10 keV released in the IE 2259-1-586 out- 
burst was 3 X 10^^ erg and 2 x 10'*° erg during the rapid 
and gradual decay components, respectively (see Table 5 
of Woods et al. 2004). For IE 1048.1-5937, we estimate 
that up to 2007 April 28, the date of the last CXO obser- 
vation, the total energy emitted is comparable, roughly 
7 X 10'^^ erg (2—10 keV), based on the pulsed fraction- 
total flux correlation and assuming a distance of 2.7 kpc 
(Gg,cnslcr et al. 20 0j), or 8 x 10 ^° erg assuming 9.0 kpc 
( Durant fc van Kerkwiikll2006af) . Near-IR flux variabil- 
ity, possibly (in the case of IE 1048.1-5937) and hkely 
(in IE 2259-1-586) correlated to X-ray flux, was observed 
following both events. That the 2007 IE 1048.1-5937 
event clearly involved a large rotational glitch, as did 
the 2002 IE 2259-1-586 event, is also a commonality, al- 
though glitches at the time of the 2001-2002 flares could 
easily have been missed because of the large amount 
of timing noise and our inability to phase connect the 
data at that time (Dib et al., in preparation). SGR- 
like bursts from IE 1048.1—5937 have not been observed 
thus far in 2007, but that does not preclude their exis- 
tence given that our observing duty cycle is so low. Sim- 
ilar to IE 1048.1-5937, IE 2259-^586 exhibited a pulse 
profile change, albeit one involving only the fundamen- 
tal and first harmonic, in constrast to the appearance of 
higher harmonic structure in IE 1048.1—5937. Finally, 
as in IE 1048.1—5937, spectral hardening was observed 
during the IE 2259-^586 flare. Thus overall, the 2007 
IE 1048.1—5937 event has practically all properties con- 
sistent with the 2002 IE 2259-F586 event, but differs from 
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its earlier flares primarily by rise time. 

The transient AXP XTE J1810-197 presents an- 
other interesting example for comparison. This 
AXP, which "turne d-on" in late 2002 or early 2003 
(jibrahim et al.l 120041 ) . has been steadily declining in X- 
ray fl ux since, reaching appare nt quiescence in 2005- 
2006 (|Gotthelf fc Halpernl I200I . Unfortunately, very 
little is known about the onset of enhanced emission, 
ie., whether it was accompanied by a glitch/burst, or 
what the timescale of the rise was. While the ex- 
ponential decay timescale of several years does resem- 
ble that of IE 1048.1-5937's 2001-2002 post-flare be- 
havior, how the most recent flare will fade is as yet 
undetermined. Furthermore, the peak and quiescent 
fluxes of XTE J1810-197 differed by a much greater 
amou nt, nearly 2 orders of magnitude (jGotthelf et al.l 
|2004[) . while the recent total 2-10 keV flux increase 
in IE 1048.1-5937 was by a factor of ~7. Spectrally, 
XTE J1810— 197 was harder in outburst and softer in qui- 
escence, similar to IE 10 48.1-59 37, when modelled as 2- 
temperature blackbody (jGotthelf fc Halpcrn 2007) or as 
a blackbody plus power law (W. Zhu, private communi- 
cation). One major distinction, however, is the clear ev- 
idence for a positive correlation bet ween pulsed fraction 
and total flux in XTE J1810-197 f Gotthelf fc Halpern' 
l2007f ). While it may be tempting to link the transience 
of XTE J1810-197 to IE 1048.1-5937-like post-flaring 
behaviour, the dissimilarities are important and remind 
us that we may be comparing mere coincidences. 

Compared to the SGR giant flares, rare events in 
which an enormo us amount of bro ad band energy 
(10'''* - 10^6 gj.g. iHurlev et all [2005h is output, the 
IE 1048.1—5937 event is orders of magnitude less en- 
ergetic. Giant flares release most of their energy in a 
short (<1 s), initial 7-ray spike, a property not observed 
in this event, although a small initial burst could have 
been missed. Furthermore, the 2007 event is prolonged, 
and while giant flares do gradually decay, the amount of 
energy in the tail is usually small compared to that of 
the spike (IWoods et al.ll2004l : llmlev et al.ll2005[ ). 

On the other hand, the gradual evolution of the persis- 
tent properties seen before and after the 2004 December 
giant flare of SGR 1806—20 do show some resemblance 
to IE 1048.1— 5937's behaviour. During this pre-flare 
period of enhanced burst activity in SGR 1806—20 that 
began in mid-2003, the pulsar torque, pulsed flux, to- 
tal flux, and hardness increased on a ~year timescale, 
peaking several months before the giant flare occurred in 
2004 Decem ber, and continuing to decline well after the 
flare epoch (jWoods et al.l 120071 ). The gradual changes, 
particularly the well-resolved factor of >2 increase in 
total unabsorbed flux, perhaps resemble the 2001-2002 
events of IE 1048.1-5937 more than the 2007 event. 
As with IE 1048.1—5937, a correlation was found be- 
tween spectral hardness and intensity in SGR 1806—20, 
and there appeared to be some sort of physical con- 
nection between frequency derivative (or torque) evo- 
lution and variability in the spectrum and pulsed flux 



(jMereghetti et al.ll2005HWoods et al.ll2007D . However, in 
SGR 1806—20, there was no evidence for a correlation 
between the phase-averaged flux, which likely peaked 
around 2004 October, and pulsed fraction, which was sta- 
ble in the flux-enhanced pre-flare period but was lower 
than average in early 2005 when the persisten t flux had 
approached quiescent levels (|Rea et al.ll2005bf ). Further- 
more, SGR 1806— 20's pul se profile was notic eably more 
sinusoidal while brightest (jWoods et al.|[2007l ). although 
the pulse profile likely changed as a direct results of the 
giant flare, rather than the slow evolution of persistent 
flux. 

4. CONCLUSIONS 

We have considered ~10 yrs of multiwavelength obser- 
vations of IE 1048.1-5937. The source, in 2004-2006, 
gradually relaxed apparently back to quiescence in most 
respects, only to "awaken" suddenly in the glitch/flare 
event of 2007 March. By observing IE 1048.1-5937 
while quiescent, it has become apparent that previous 
instabilities in timing. X-ray and near-IR flux were likely 
all linked to the major long-term flaring events of 2001- 
2002, whose nature is yet unknown but likely includes 
changes in the stellar magnetosphere's current configu- 
ration. The asynchronous spin and X-ray flux variability 
we have observed is incompatible with expectations of 
the fossil disk accretion model. Following the most re- 
cent 2007 event, we observed total X-ray flux variability 
that is strongly correlated with X-ray pulsed fraction, X- 
ray spectral hardness, and changes in near-IR flux. These 
observations largely agree with expectations of the mag- 
netar model. To date, the IE 1048.1-5937 2007 event is 
ongoing, as is continued multiwavelength monitoring. 
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